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We argue that (first-order) coherence is a relative, and not an absolute, property. It is shown how 
feedforward or feedback can be employed to make two (or more) lasers relatively coherent. We also 
show that after the relative coherence is established, the two lasers will stay relatively coherent for 
some time even if the feedforward or feedback loop has been turned off, enabling, e.g., demonstration 
of unconditional quantum teleportation using lasers. 
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PACS numbers: 42.50.Ar, 42.50.Lc, 42.55.Ah 



I. INTRODUCTION 

Recently, several authors have pointed out that al- 
though it has been known for a long time that a laser 
does not generate a coherent state, it is still possible to 
use laser light to observe phenomena that are theoreti- 
cally described in terms of optical (first order) coherence 

Qilaa _ 

M0lmer 111 argues that optical coherence, manifested, 
e.g., in a non- vanishing expectation value of the electric 
field operator, is a "convenient fiction" and he demon- 
strates that experiments that are usually interpreted in 
terms of optical coherence need not be described in those 
terms. Moreover, he argues that optical coherence is not 
easily generated, and in particular, that a laser does not 
generate coherent state. Rudolph and Sanders '5| show 
that if M0lmer's assertions are true, so that the output 
state of a laser is given by the density operator 
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then quantum teleportation is not possible under the con- 
ditionspostulated for "unconditional quantum teleporta- 
tion" " ^ . Rudolph and Sanders do not rule out the ex- 
istence of sources of coherent states, but argue that lasers 
do not produce light with (first order) optical coherence. 
Van Enk and Fuchs Q claim that while a laser does not 
produce optical coherence in M0lmer's sense, subsequent 
temporal modes of the laser output are phase correlated. 
That is, if the output from a laser is expanded in subse- 
quent orthogonal temporal modes, the state of M such 
modes is 
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where the symbol denotes a M-fold tensor product. 
A requirement for this is that the laser's coherence time 
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is longer than M times the time duration of a tempo- 
ral mode. In such a state all the M subsequent modes 
are first order coherent relative to each other, in contrast 
to the M-mode state /5f*^. Both states have a vanish- 
ing expectation value of the electric field operator for 
all modes. Van Enk and Fuchs claim that such "phase 
coherence" is sufiicient to allow unconditional quantum 
teleportation. Wiseman Q , finally, asserts that there are 
"no devices that can generate 'true coherence' any better 
than a laser." The basis for his claim is that any oscilla- 
tor, or oscillation, can only be described relative to some 
accepted "clock" standard. He goes on stating that noth- 
ing suggests that there exist any better "clock" at optical 
frequencies than a laser. We subscribe to this view, and it 
is compounded by the fact that at the National Institute 
of Standards i Boulder, Colorado, the next generation 
of an "atomic clock" in development is indeed an opti- 
cal clock 7]. That is, the clock oscillator operates at an 
optical-, rather than a microwave-transition of an atom. 

The coherent state |e*'^|Q:|) of an oscillator with the an- 
gular center frequency loq, can mathematically be gener- 
ated from the vacuum state through the displacement op- 
erator. It is a minimum uncertainty state in the in-phase 
operator ai = (a + a^)/2 and the quadrature-phase op- 
erator 02 = (a — a^)/2i, where a (at) is the bosonic anni- 
hilation (creation) operator. The quadrature amplitude 
fiuctuation operators are defined Aa^ = di — (a^), with 
i = 1,2. In Table n we have computed the expectation 
values of the quadrature amplitudes, the photon number, 
and their variances for the two states pL and |e*'*|a|) (for 
the particular choice cf) — Q) in columns two and three. 
The table clearly shows that the observable statistics of 
the two states differ significantly. Hence, they cannot 
both describe the output of a laser, at least not if simi- 
lar initial and boundary conditions for the laser and the 
detector are assumed. Wang et al. have also studied and 
quantified this difference |1 ij . 

In the following we shall describe a simple model of 
phase locking of two independent lasers to a master laser. 
Phase-locking of two lasers has already been demon- 
strated using both CW 8] and pulsed lasers 0|. How- 
ever, phase-locking of two independent lasers close to the 
quantum limit, that is such that the relative phase fiuc- 
tuations are comparable to those of two ideal coherent 
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states, has never been demonstrated. We shall see, that 
it is possible to stabilize two lasers so that they are rel- 
atively first-order coherent to that extent, although they 
are incoherent relative to any third, auxiliary laser, or 
other "clock standard" . A similar proposal has also been 
put forth by Fujii 12], but in ^2], the analysis is geared 
towards a practical implementation of quantum telepor- 
tation, whereas we direct our interest towards the ques- 
tion whether first order coherence between two lasers can 
be established at all. 

It will be convenient to work in the Heisenberg picture 
in the spectral domain. Since a CW laser output consists 
of a continuous photon flux, we shall work with the pho- 
ton flux operator f. (See subsection IA II in the Appendix 
for a detailed discussion of f.) For a coherent state, the 
spectral relations corresponding to equations three and 
four in the third column of Table HI are ^3 

SAf^m = 1/4, (3) 
SAf^m - 1/4, (4) 

where £7 = a; — cjo so that, e.g., SAfii^) is the double- 
sided power spectrum of Afi around the center frequency 
luq. (Note that the laser external field spectra com- 
puted in Ref. [Kj are single-sided spectra.) These re- 
lations demonstrate a particular feature of the coherent 
state, namely a frequency independent quadrature am- 
plitude noise spectrum. Operationally, this means that 
the quadrature amplitude noise of a coherent state is sta- 
tionary. Moreover, a field in a coherent state remains in 
a coherent state independent of the detector temporal 
response function (see the Appendix). 

The corresponding spectra for the laser external field 
have also been computed in |0. While the model em- 
ployed was specifically targeting a semiconductor laser, 
the general features are largely independent of the laser 
type. It was found that the corresponding relations for 
(a somewhat idealized) laser, pumped high above the 
threshold, are 

SAf^{^) = 1/4: and (5) 

Sa,A^^) = (6) 

where 7 is the inverse of the laser photon lifetime. We see 
that © radically differs from Q in that Eq. lO has a 
l/il^ behavior at low frequencies. This is characteristic of 
a Wiener process, leading to diffusion. Therefore, even 
if the phase of the external field amplitude was known 
(e.g., through a series of measurements) at some time to, 
the phase will be randomized through phase diffusion at 
times t ^ to + I/7. Therefore, unless the field of the 
laser has been measured relative to a reference within 
the laser's coherence time I/7, the state of a laser is well 
described by the density operator irrespective of how 
relatively stable the reference oscillator is. The density 
matrix also describes the state of a free-running laser, 
that is, a laser whose phase relative to our reference is 
unknown. This result is indeed expected from symmetry 



arguments, and M0lmer, in particular, emphasize that 
there is no mechanism in a free-running laser that breaks 
the in-phase, quadrature-phase symmetry. Hence, the 
expectation value of the electric field of an unmonitored, 
free-running laser vanishes. 




FIG. 1: A schematic drawing of two lasers locked by feedfor- 
ward to a master oscillator laser. 



II. PHASE-LOCKED LASERS 

Now consider the case when we stabilize the laser field 
with the help of a master oscillator. In the following 
we will consider feedforward stabilization, for simplicity. 
Later, we shall briefly discuss feedback stabilization. Fig. 
n shows a schematic drawing of two lasers locked to a 
master oscillator by feedforward. The notation we will 
use is defined by the figure. The master oscillator (laser) 
field b is divided into two by beam splitter one (BSl) 
which is a 50/50 beam splitter. The field into the other 
port of the beam splitter is a vacuum field u, so that 
the beam splitter output modes k and k are given by 
k = -^=[b + u) and k = ^ • 

Next, we look at beam splitter two (BS2), where part of 
the laser field a is tapped to get a locking signal d. Again, 
a vacuum field, denoted v, is incident on the other input 
port. Hence, the two output fields d and c are given by 

d — rd + ti), and c = td ~ rv, (7) 

where r = sin 9 (t = cos 9) is the reflectivity (transmis- 
sivity) of the beam splitter. At beam splitter number 
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three (BS3), which is also a 50/50 sphtter, the fields d 
and k are mixed, and subsequently the outputs g and h, 
which are given by 
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^ {d + k)/V2^ {ra + tv)/V2 + l{b + u) (8) 



and 



{d~k)/V2= {rd + tv)/V2-^{b + u) (9) 

are measured by photodetectors. The detector photocur- 
rents will be proportional to g^'g and h^'h, that can be 
expressed as 
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^g= {d^d + Sk + k''d + k^k)/2 



and 



h^h = {d}d-d)k - yd + k^k)/2. 



(10) 



(11) 



The difference between the measured photocurrents g'^ g~ 
h^h will be the homodyne signal e, which we call the error 
signal. This signal is simply 



= Sk + k''d. 



(12) 



Now, expand the operators in quadrature components 
and assume that d leads k by the relative phase angle 
7t/2. (this choice of relative phase will beat the excita- 
tion of one field with the quadrature amplitude of the 
other field.) That is, we multiply the operator d with 
exp(i7r/2) = i and remember that both fields are ex- 
pressed relative to the same fiducial field. (Below, we 
shall see that the relative phase angle 7r/2 is the stable 
feedforward locking point.) This will yield the mean error 
signal zero, and the fluctuating part of the error signal 
will be 

e = -i{do + Adi - iAd2)iko + Aki + iAh) (13) 
+i{ko + Aki - iAk2)ido + Adi + iAd2) (14) 
= 2{doAk2-koAd2) + 0{A^), (15) 

where ©(A^) denotes terms of second order in the 
quadrature flux fluctuation operators. 

It is well worth examining Eq. H15|) in some detail, 
since it is the generic noise equation for homodyne mea- 
surements. Assume, that each of the fields k and d 
is in a coherent state. Hence, it follows that the two 
fields are relatively coherent. Since the quadrature am- 
plitude fluctuations of a coherent state are equal, and 
independent of the state's excitation, the homodyne de- 
tection noise is dominated by the weaker field, since the 
weaker field's quadrature amplitude noise beats against 
the stronger field's in-phase amplitude. (It is custom- 
ary to refer to the stronger field as the "local oscilla- 
tor.") If the two fields have unequal quadrature ampli- 
tude fluctuations, e.g., assume that S^f.^{n) > S^^^{Q), 
then the homodyne detection noise is still dominated by 



the quadrature amplitude fluctuations of d provided that 

kl>dls^f^^{n)/s^^^m. 

Using the relations between incident and output fields 
above, we can re-express this equation in the (quadrature 
expansion of the) input fields a, b, u, and v. The result 
(to first order in the fiuctuation operators) is 
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This is almost what we desire. Remember that the ob- 
jective is to reduce the large fluctuations of A 02 at low 
frequencies. These are the fluctuations leading to phase 
diffusion. The error signal clearly contain the information 
about these fluctuations. By making the master oscilla- 
tor flux amplitude sufficiently large, so that 60 is much 
larger than oq, this noise information will dominate the 
error signal at low frequencies. (Remember that the spec- 
tral density of Aa2 and A62 increases as whereas 
the spectral density of Au2 and Av2 is 1/4, indepen- 
dent of frequency. ) From the expression IjlGI) we see that 
under this condition, if the relative phase between Aa2 
and A62 is larger than tt/2, that is, either Ad2 is posi- 
tive or A62 is negative, then the error signal is negative. 
In order to compensate for the Aa2 fluctuations, the er- 
ror signal should be feed forward, keeping the sign, to a 
phase shifter. Let us, therefore, see what happens when 
we apply some phase shift, say i^, to an output signal 
c = Co + Aci + iAc2. In this case we denote the phase 
shifted signal by / 

/ = exp(i(/))c (17) 

= (co -I- Aci) COS0 — Ac2 sin0 (18) 

+i[Ac2C0S(j)+ {co + Aci)sm(f>] (19) 

« co + Aci+i(Ac2 + co0), (20) 

where we have linearized the equation to first order in 
(f>, Aci, and Ac2. Assume now that the phase shift (/) is 
equal to the error signal e times some feedforward gain 
Q. The output signal after phase correction by the error 
signal is given by 

/ w t{ao + Aai + iAd2) - r{Avi + iAv2) + iiaot/e.(21) 

To cancel the noise term Aa2 in the expression above, 
we have to adjust the feedforward gain so that 



G = V2/raobo. 
After doing so, the output signal becomes 

tao , 



(22) 



/ = tao + t Afli — rAvi + i 



-(A62 + AU2) 



-AV2 



(23) 

where we have used the relation P + r'^ = 1. From this 
equation we see that the noise term Ad2, describing the 
phase-diffusion noise originating in the laser generating 
the d field, is absent. Instead, a new diffusion noise term 
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taoA62/&o has appeared. Since the master oscillator is 
also a laser, this noise term will also have a spectral den- 
sity proportional to il^^ at low frequencies. Hence, one 
Wiener process, Aa2 (relative to the fiducial reference), 
has simply been replaced by another, A62- 

Let us now see what happens with the output field / 
if we subject the laser field a to the same sequence of 
measurement and feedforward as the field a. The calcu- 
lation of / proceeds in a similar fashion as for /, except 
that the field k = (h — u)/\/2. Hence, the only difference 
from the expression H23|) above is that all fields should be 
primed, except for the master oscillator field h and the 
vacuum field operator u. However, the sign in front of 
the operator u should be reversed. Hence, the output / 
becomes 



of that of a coherent state. For a more detailed discus- 
sion of the connection between (|25|l and the first order 
coherence function between the locked lasers, the reader 
is referred to the Appendix, subsection I A 31 

An additional consequence, discussed in more detail in 
subsection lA 21 in the Appendix, is that within a time 
small compared to the inverse spectral linewidth of the 
locked lasers, the lasers will remain relatively coherent 
even if the feedforward locking loop is turned off. It will 
thus be possible to do, e.g., continuous variable teleporta- 
tion without giving up the quantum teleportation criteria 
given in 0| 



III. FEEDFORWARD V.S. FEEDBACK 



bo r 



(24) 

Assume that the two field amplitudes tao and tag are 
equal and that the two fields / and / are detected by a 
balanced homodyne detector. The detector output, that 
is the beat-note between the two incident fields / and 
/ , can then be deduced in a straightforward fashion by 
comparing with H15|l . The detector output fluctuations 
will be 



/oA/2-/oA/2 = 2tao 



2taQAu2 Auj — Ail 



(25) 



From this equation we see that the relative-noise spec- 
trum is no longer proportional to fl^^ at low frequencies, 
but is flat, since Au2, Av2 and Av2 all emanate from vac- 
uum fluctuations and hence have frequency independent 
spectra with the spectral density 1/4. Recalling that an 
earlier assumption was that bo 3> ao, the first term on 
the right hand side of (|25|l can be neglected. Therefore, 
the spectral density of the detector fluctuations is, to 
a good approximation, 2(iao/r)^ at all frequencies (see 
the Appendix). In comparison, the homodyne detector 
spectra from two coherent states, with field amplitudes 
too, would be 2(too)^. Hence, as r ^ 1, the relative 
quadrature-phase noise of the two locked lasers will ap- 
proach the noise level of two coherent state sources. At 
the same time, the mean in-phase amplitude of / and 
/ would approach zero. Assume now that we use the 
relatively coherent fields / and / as local oscillators in 
two separate homodyne measurements. Using (I15|l . it is 
not difficult to prove that in order to minimize the influ- 
ence of the local oscillator quadrature amplitude noise, 
the choice t = r = 2~^/^ should be made. Assume, for 
simplicity, that the field we want to measure is a coherent 
state with mean amplitude do. If so, the coherent state's 
quadrature amplitude noise will dominate the measure- 
ment fluctuations if ao > 2do. That is, with this choice 
of mirror transmission t, the two flelds / and / are rela- 
tively coherent, each with a noise spectral density twice 



Above, we have analyzed the situation where two lasers 
are phase locked to a third laser by feedforward. While 
such a scheme lends itself to a simple analysis, the scheme 
has obvious practical shortcomings. One is the need of a 
precise control of the feedforward gain. In order to com- 
pletely suppress the phase diffusion of each laser, relative 
to the master laser, the feedforward gain must be pre- 
cisely set according to (|22ll . This will be impossible in 
reality. In addition, the mean absolute value of the error 
signal (|I6|I will be grow with time as t^^^, there t is the 
time since the main information the error signal contains, 
Aa2, is a non-stationary term. The error signal is used to 
drive a phase-compensating device, e.g., an electroopti- 
cal crystal or an adjustable piezoelectrically controlled 
delay line. However, due to the fact that with time, 
the error signal will increase without limit, any phase- 
compensating device will eventually run out of range and 
will have to be reset. 

In reality, it would be wiser to lock the lasers by feed- 
back, and the feedback should act on the laser frequency 
(e.g., translating one of the laser's mirrors) after appro- 
priate filtering. In a feedback loop, the relative fiuc- 
tuations will initially become smaller and smaller with 
increasing feedback gain. If the feedback measurement 
and feedback loop is carefully designed, it is possible to 
reach the limit set by quantum noise manifested in H25() 
before (non-fundamental) feedback loop fluctuations are 
sufficiently amplified to dominate the locked lasers' out- 
put. In addition, since the accumulated relative phase 
equals the time integrated frequency difference between 
the master and slave laser, the phase-compensating de- 
vice will not run out of range. Differently stated, the 
lasers frequency difference is the time derivative of the 
lasers' relative phase. Therefore, if the latter noise pro- 
cess has a spectrum proportional to i^^^, the former noise 
has a fiat spectrum at low frequencies. The frequency er- 
ror signal is hence a stationary process. Hence, if the 
feedback is implemented by moving one of the laser's 
mirrors, its position needs only to be adjusted within a 
fixed range and the actuator moving the mirror need not 
run out of range. 

The disadvantage with feedback locking is that it is 
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more difficult to analyze, and that one risks self oscilla- 
tion at baseband frequencies Q where Orj = it. Here, 
Tf is the feedback loop time delay. At this frequency the 
feedback loop does no longer suppress the fluctuation but 
instead enhances them. If the feedback loop gain is too 
high, one induces self oscillation. Hence, one needs to 
make a compromise between the feedback loop stability 
and the feedback gain, all while maximizing the feedback 
loop bandwidth. Since our analysis essentially only had 
the purpose to point out that coherence of a harmonic 
oscillator must always be relative to some reference os- 
cillator, we will not delve deeper into these issues here. 



IV. CONCLUSIONS 

We have argued that first order coherence is a relative 
property, and that if one accepts this view, the light em- 
anating from a laser can be first-order coherent relative 
to itself at an earlier time, or to a different laser. Due 
to phase diffusion, two lasers will only stay relatively co- 
herent for a time smaller than the smallest of inverses 
the lasers' respective linewidth unless they are actively 
phase-locked. Since, in principle, the laser linewidth can 
be as small as one wishes, this fact does not prevent lasers 
to be used to demonstrate, e.g., unconditional quantum 
teleportation. It could in principle be done by locking 
Alice's and Bob's homodyning lasers to a master laser 
before the unknown state and the shared entangled are 
measured with the help of the laser. The locking is then 
turned off, and Ahce's and Bob's respective laser become 
free running. The two experimentalists now have roughly 
the lasers' coherence time to perform the teleportation, 
including the Alice's homodyne measurements, the clas- 
sical communication, and the final unitary evolution by 
Bob. A typical laser's coherence time may seem short on 
the human time-scale (~1 s). However, on an electron- 
ics time-scale (~1 ns), a two good lasers stay relatively 
coherent for a long time. 

In the Appendix, we have included a discussion why 
the popular model of laser fluctuations using Langevin 
noise sources fail to describe the laser's behavior for 
times long compared to the laser's inverse linewidth. 
The linearization customarily employed in the solution 
of these equations neglect to take the restoring force of 
the quadrature-phase amplitude fluctuations in account, 
and hence these equations erroneously predicts that it is 
only the quadrature-phase amplitude, relative to some 
fix fiducial reference, that diffuses. Unfortunately, if the 
equations are not linearized, they are difficult to solve 
analytically. 
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APPENDIX A: 



EXPECTATION VALUES FROM 
FLUX SPECTRA 



1. Expectation values of a coherent state 

Here we derive the relation between the external field 
operator f and the mean photon number and the photon 
number variance of a specific mode. The external field 
operator f correspond to a photon flux amplitude. The 
photon flux operator is hence given by f^f, and has the 
unit Hz. The expectation value of f^f gives the mean 
photon flux, in photons per second. Let us expand f 
relative some fiducial reference as f{t) ~ vq + Afi(f) + 
i Af2(t). Here, the field amplitude tq has been assumed to 
be in phase with a fiducial reference. Expanding to first 
order in the fluctuation operators Afi(i) and Af2(t) we 
get 



(Al) 



Assume now that a photodetector temporal mode- 
function /i is a rectangle function of duration T, see Fig. 
121 The detector output, as a function of time, will be the 
convolution between the flux and the detector response 
function 



R{t) = {f^f) *hx[r^ + 2roAfi(i)] * h. 



(A2) 



hit)k 



1 













-T/2 



T/2 



FIG. 2: The assumed temporal mode function. 

Sampling this function every T seconds, yields the pho- 
ton number in subsequent orthogonal temporal modes. 
The detected mean photon number per mode hence be- 
comes 



(A3) 



We would now like to compute the photon number vari- 
ance of such a mode. Since the coherent held flux fluc- 
tuations are ergodic we can proceed as follows: Define 
y{t) = Afi(i) * h. The spectral density of this signal is 

= Sf,it)\Hm\' - ]\Hm', (A4) 
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where the spectral filter function H{^) is the Fourier 
transform of h and we have used the fact that the spec- 
tral density of the quadrature amplitude fluctuations of 
a coherent state flux is 1/4 according to (|2Jl. The first 
order correlation function G{t) of y(t) is the (inverse) 
Fourier transform of its power spectrum, so that 



-h * h, 



(A5) 



where we have used the fact that the Fourier transform 
of a product of functions is the convolution of their re- 
spective inverses. Finally, the variance ofy{t) is Gy(t)(0). 
Looking at figure |3 we find that Gy(t) (0) = T/A. This 
means that the photon number variance can be expressed 

{Afi^) = 4rl{{Ah{t)) * hf) = ArlT/A = Tr'l (A6) 



Comparing with (|A3p we see that we arrive at the ex- 
pected result for a coherent state. 




FIG. 3: The temporal mode function convolution. 

To derive the mean quadrature amplitudes and their 
variances we proceed in a similar way. Suppose two fields, 
characterized by the photon flux amplitudes f and f are 
incident on a homodyne detector. If the two fields are 
relatively coherent and have the same phase, we can use 
(|12|l and insert the expansion of the fields in quadrature 
amplitude components (e.g., f = + Afi + iAf2). To 
the first order in the fluctuation operators the detector 
output becomes: 



2roro + roAf^ + r^Ah + OiA^). 



(A7) 



Passing this signal trough appropriate detector filter, (in- 
tegrate for T seconds and then dump), the mean detec- 
tion signal becomes 



2rorgT. 



(A8) 



The prefactor 2 in the equation above is simply a scaling 
factor. Since the signal emanates from the projection of 
one in-phase amplitude on another, we conclude from the 
symmetry of the problem that, for the field f. 



(A9) 



(Remember that, e.g., ai is the in-phase amplitude of the 
field in a specific temporal mode, while fi is the in-phase 
amplitude operator of a continuous flux.) Assume now 
that f has a much larger mean in-phase amplitude than 
f. Filtering away the DC-component of the signal l|A7|l . 
and measuring the variance of the detector signal after 
flltering the temporal mode, the detector signal approx- 
imately becomes 



TgAfi * h{t) 



(AlO) 



Using the same procedure as we used to compute the 
mean photon flux above, we see that the variance of 
this signal is {rQ)'^T/4:. From this we can conclude that 
{AdD = 1/4, since T{rQ)'^ is the square of the mean local 
oscillator in-phase amplitude (the amplitude we beat the 
fluctuations against). 

If we phase rotate f with tt/2 with the respect to f , and 
follow the same procedure to calculate (0,2) and (Aaf), 
we will find that (0,2) = and {Ad\) = 1/4. Hence, the 
results in the third column of Table are reproduced if 
Eq. (O and igj) are used. 



2. Expectation values of the laser output field 

If we now use (jSJ and modelling the noise prop- 
erties of the output field from a laser, we will get iden- 
tical results as those in the previous subsection for (n). 



{{Ahf 



{ai) 



{a2) 



and (Ac 



However, due to the 



spectrum Q of f2, the expectation value of (Aa^) = 00 
independent of T. 

To study this divergence, due to the diffusion of f2, 
assume that we homodyne the field f in two orthogonal 
temporal modes defined by two rectangle functions with 
durations from t to t+T, and t+r to t+r+T, respectively, 
where t > T. A schematic procedure to make such a 
measurement with the aid of a flip mirror is depicted in 
Fig. ^ From (|15|l we find that the measurement output 
wiU become 2ro[Af2(t) - Af2(i + t)] * h{t). With the 
help of ©, the spectrum of the output fluctuations can 
be computed to be 



sin(17r/2)sin(r2T/2) 



n/2 



(All) 



Transforming this function back to the time domain gives 
the correlation function, and its value at the origin, the 
variance, can be computed to be: 



2-rT7^(i-f ) 

or 



(A12) 



We see that since t >T, the variance will approximately 
be ^TqT (the shot noise associated with the total detected 
photon number ir^T) when Tt < 7"^. That is, for mea- 
surements, even relative measurements, within the laser 
coherence time 7""'^, the laser output held in a succession 
of temporal modes will, each one, approximately be in a 
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coherent state if the first mode in the succession is used 
as the reference. Note that this conclusion holds irre- 
spective of the numerical values of T, r, and 7"^ as long 
a,s T < T < . This result agrees with the analysis of 
van Enk and Fuchs 0. 



3. First order coherence of phase-locked lasers 

The first order correlation function between two clas- 
sical electromagnetic fields is defined 



{E{f,t)E'*{f',t')), 



(A13) 



where r and r' denote two spatial locations, and t and t' 
denote two times. In the following we will restrict our- 
selves to single mode fields, and therefore, we can sup- 
press the spatial coordinates. Expressing the fields in 
terms of amplitude and phase, and expanding the ampli- 
tude in its (real and positive) mean Eq, and its fluctua- 
tion around this mean AE{t), one gets 

tujo{t-t')/i[4>{t}-4,'it')]\ 



{E{t)E'*it'))^EoE^e 



(A14) 



where we have dropped the term quadratic in the the 
fluctuations AE, and where we have assumed that the 
amplitude- and the phase-fluctuations are uncorrelated. 
(The latter assumption is not true in semiconductor 
lasers, where the so-called a-parameter characterizes the 
inverted-media meditated correlations between ampli- 
tude and phase. For small fluctuations this fact is in- 
consequential for what follows.) From Eq. (|A14|I we see 
that as long as \(j>{t) — (p'i^')] ^ 1 rad, then the corre- 
lation function has the modulus \EoE'q\, indicating that 
the two fields are first order coherent. Assume that the 
diff'erence \(j){t) — (/>'(t')l is small. We can then take either 
field and make it our fiducial reference, so that we work 
in a frame rotating with the angular frequency luq- We 
find that in this rotating frame, Eq ^ tq, AE{t) ^ Afi, 
and for small angles we have (f)(t) <-> Af2/ro. Similar re- 
lations hold for the primed field. Hence, if we can show 
that for two locked lasers the relation 



\Af2/ro~Af'^/r'o\ « 1 



(A15) 



holds, then our assumptions above hold, and the laser 
fields will be first order coherent. For two lasers with 
equal fields, rg = Tq, the condition can be reformulated 

|Af2- Af^l <ro. (A16) 

Hence, the homodyne measurement beat signal is directly 
relevant to the two fields' relative first order coherence. 

From H25|l we can deduce that the noise spectrum of 
Af2 — Af2 equals 4, where we have assumed a mirror 
reflectivity of 1/2 for the beam splitters BS2 and BS2'. 
This corresponds to an equivalent photon flux amplitude 
of the order 1 Hz^/^, very much below a typical laser. 
E.g., a single (transverse) mode laser emitting 1 mW op- 
tical power at a wavelength of 500 nm has a photon flux 
amplitude of 5 • 10^ Hz^/^ demonstrating that ljA16(l is 
easily met. This proves the validity of our approach. 



Flip mirror 

'%, 

iLaserT >^»<(^- 




-T/2 T/2 



FIG. 4: Schematic drawing of a homodyne measurement of 
different temporal modes of a laser's output flux. The flip 
mirror is synchronized with the detector, and the separation 
between the temporal mode can be varied through the time 
delay t = 21 /c, where I is the length of the time-delay "trom- 
bone" and c is the phase-velocity. 



4. Laser potentials and the linearized 
approximation 

To see why the noise spectra , ^ correctly predicts 
the short-time (T, r < 7"^) relative-coherence properties 
of a laser, but fails to reproduce the long term quadrature 
phase properties (popularly speaking, these are called the 
coherence properties) manifested in the flrst column in 
Table it is instructive to look at the laser potential 
models that underlies the various theories. 

The internal field of a laser, under certain approxi- 
mations can be modelled as a particle in a particular 
potential subjected to Langevin noise forces. The sim- 
plest standard model for the potential, derived from the 
equation of motion of the field is 



(A17) 



where a is the overall gain per unit time in the laser, 70 
the cavity decay rate (the inverse of the (cold) laser cavity 
photon lifetime), and C is the gain saturation parameter. 
In Fig. |31we see how the potential goes from a parabolic 
potential to a "Mexican hat" potential as the gain goes 
from subthreshold to above threshold. In both cases the 
potential has circular symmetry around the origin. Since 
lasing is triggered by spontaneous emission (that is one of 
the origins of the Langevin noise sources in the model), 
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FIG. 5: A laser potential with 70 = 1 per unit time and 
C = 0.8. (a) Below threshold pumping {a — 0). (b) Above 
threshold pumping (0 — 70 = 1.6 per unit time). 

the laser field will not have any preferred phase (rela- 
tive to the fiducial reference defining the coordinate sys- 
tem orientation). Therefore, a free-running laser where 
only the intensity is known (supposedly through a pho- 
ton number measurement or by knowledge of the relation 
between the pumping and the output intensity) must be 
described by the density operator /5l of Eq. QJ. As 
M0lmer P|, and other's before him, have pointed out, a 
laser does not induce any symmetry breaking, and there- 
fore it does not induce coherence relative to any other 
oscillator. 

However, as the analysis above indicate, by measuring 
and influencing the field of a laser, either with feedfor- 



ward control, or feedback control, the relative phase of 
two lasers can be locked so that the lasers become rela- 
tively coherent. However, if the locking servo is turned 
off, the Langevin fluctuations of each laser will make the 
relative phase between the two lasers diffuse, since the 
potential offers only a restoring force in the radial direc- 
tion, but not in the azimuthal direction. After a time 
approximately equal to the inverse of the laser emission 
spectral linewidth, the two lasers are no longer relatively 
coherent. 

This situation is not well described by Eqs. lO and 
© , because these equations predict that it is only the 02 
quadrature that diffuses. If that were the case, then the 
relative phase of the two lasers with the initial relative 
phase zero, would remain zero, on average, for all subse- 
quent times. If so, the quadrature-phase amplitude fluc- 
tuation would eventually become larger than the mean 
in-phase component. However, this is not what happens 
in a laser. The reason for the erroneous prediction of 
© is that in the analysis leading to the equation, the 
laser potential (or rather the Langevin noise operators) 
have been linearized about the operating point. The cor- 
responding linearized potential is illustrated in Fig. |S1 
One can see that in the linearized potential, the gain sat- 
uration will only act in the direction parallel to the field 
expectation value. In contrast, the better model of the 
potential will act in the direction parallel to the instan- 
taneous field, i.e., in the radial direction. The linearized 
potential correctly predicts the laser behavior for times 
up to the laser's relative decoherence time, but fails for 
long time scales, where it is actually the relative phase, 
and not the quadrature field amplitude that diffuses. 




FIG. 6; A linearized laser potential above threshold 
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PL 16**^0) Laser model 

(ai) |q| ^/Tro 

(aa) 

(Aa?) ^ii£±i 1/4 1/4 

(Aai) ^ 1/4 

{n) = (ata) |Qp jaj^ Trg 

{(An)2) lap lap Trg 



TABLE I; A table comparing the predicted statistics of three 
different fields. A totally statistically mixed state with Pois- 
sonian statistics, a coherent state, and the predicted output 
of a laser using Langevin equations. 
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